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ABSTRACT 

Electronic  contributions  to  the  optical  ( hyper )polarizabilities  of  small  silicon  clus¬ 
ters  are  theoretically  determined.  Geometries  and  the  electronic  structures  of  the  clusters 
are  established  using  the  tight-binding  model.  The  nonlinear  polarizabilities  are  found 
to  depend  primarily  on  the  symmetry  of  the  cluster  and  prove  to  be  high  for  the  low- 
symmetry  clusters.  Possible  experiments  and  applications  Eire  discussed. 

INTRODUCTION 

Small  semiconductor  clusters  in  the  range  from  a  few  atoms  to  tens  of  atoms  are  of 
a  great  interest  from  the  viewpoint  of  both  fundamental  science1-5  and  applications.5-8 
Their  physical  properties  (symmetry,  electronic  structure,  optical  spectra  and  transition 
probabilities)  differ  significantly  from  those  of  the  solid  state  both  in  the  bulk  and  at 
the  surface,  and  also  from  the  properties  of  nanoscale  structures,  such  as  quantum  dots. 

Linear  and  nonlinear  optical  properties  of  the  latter  have  been  shown  to  depend  strongly 
on  their  size  in  the  region  of  quantum  confinement  (see.  e.g.,  Refs.  5.  9  and  10).  But 
in  these  and  similar  works,  the  bulk  electronic  structure  of  the  semiconductor  is  usually 
assumed  and  modelled  by  free  electrons  with  effective  mass.  Obviously,  this  approach 
is  valid  only  for  sufficiently  large  objects  with  sizes  not  less  than  a  few  nanometers, 
containing  on  the  order  of  1000  atoms  or  more. 

In  the  present  work,  the  optical  properties  of  small  silicon  clusters  with  7  —  13 
atoms  are  predicted.  For  such  clusters,  the  bulk  approximation  is  not  valid  and  the 
detailed  structure  becomes  important.  Much  work  has  been  done  on  the  structure  of 
silicon  clusters,  both  experimentally  for  stability  and  photofragmentation1’2  and  for  op¬ 
tical  absorption,3,4  and  theoretically  with  ab-initio- type  calculations  for  smaller11-13  and 
other  methods  for  larger14-19  systems.  The  approach  we  use  in  this  paper  is  based  upon 

the  semiempirical  tight-binding  (TB)  model.  This  model  has  previously  been  used17  to  '  _ 

describe  the  structure  of  Siio  isomers.  It  is  relatively  simple,  thus  allowing  the  global 
optimization  of  the  geometry  even  for  comparatively  large  clusters.  Moreover,  the  TB  □ 

model  couples  the  geometry  to  the  electronic  structure,  which  we  find  essential,  but  which  □ 

usually  is  ignored  for  larger  systems.1415.  ' - 

With  the  TB  geometry  and  electronic  structure  and  using  one-electron  density 

matrix  techniques,  we  obtain  closed  sum-over-one-electron-states  expressions  from  which  - 

linear  and  nonlinear  optical  polarizabilities  of  clusters  axe  subsequently  computed.  These  _ 

characteristics  govern  a  number  of  observable  effects:  light  scattering  and  absorption  by  y  Codes 
clusters,  second-harmonic  generation  (SHG),  optical  rectification,  birefringence  induced  and/or-- 
by  optical  fields  and  the  Kerr  effect,  phase  conjugation,  etc.  lal 

THEORY 

We  write  the  TB  Hamiltonian  as18,20 


(1) 


RtB  ^  ~a  _r  ^  ^  fia  ub  0.palll,b 

Ha  pa  ub 

where  ajia  and  apa  are  electron  creation  and  annihilation  operators  in  the  basis  j fxa)  = 
^a(r  —  Rp),  with  <^>a  =  {3s.  3pz ,  3py,  3p: }  as  the  valence  orbitals  of  silicon  atom  at  the 
sites  /i  with  coordinates  Rp.  To  reproduce  the  bulk  silicon  band  structure  with  the 
nearest-neighbor  distance  2.35  A.  the  diagonal  matrix  elements  of  the  Hamiltonian  (1) 
were  fitted20  to  the  values  e°  =  -5.25  eV  and  =  1.20  eV.  and  the  off-diagonal  elements 
to  V3StT  =  -1.938  eV,  Vsper  =  1.745  eV.  \’ppa  =  3.050  eV  and  Vpp7,  =  -1.075  eV.  The  off- 
diagonal  matrix  elements  were  taken  to  behave  Slater-Ivoster-like21  in  their  angle  and 
distance  dependence,  and  so  their  bond  length  dependence  is  1/r2  up  to  3.3  A.  where  we 
consider  the  bond  to  be  broken. 17 

The  diagonalization  of  the  Hamiltonian  gives  the  one-electron  energies  sp  and  the 
eigenvectors  Cp.  Thus,  we  can  write  the  one-electron  states  of  the  occupied  and  unoc¬ 
cupied  valence  levels  as 

w  =  •  (2) 

and  the  cohesion  energy  of  the  X  atoms  due  to  the  bond  formation,  which  we  call  the 
"band  structure"  energy,  can  be  written  as 

EBs  =  Y.vp-xY.n°-°  +  uT.fi'‘-'il'>2  ■  131 

P  a  p 

where  np  and  n°  are  occupation  numbers,  and  the  difference  between  the  two  first  terms 
is  the  stabilization  energy  of  the  four  free-atom  valence  levels  due  to  the  bonding.  The 
third  term  in  Eq.  (3)  is  added  to  include  the  intra-atomic  Coulomb  repulsion  caused  by 
charge  transfer  within  the  cluster,  and  it  is  evaluated  using  the  Mulliken  charges  qp  and 
<y°,  with  the  constant  U  taken  to  be  1  eV.18 

A  repulsion  energy  term  is  finally  added  to  account  for  the  structure  of  small  clus¬ 
ters.  It  is  the  sum  of  interatomic  pair  potentials  Ed(RPv)  and  a  term  depending  on  the 
bond  number  Xb, 


.v 

P<v 

This  term  has  been  fitted18  within  the  present  TB  model  to  reproduce  the  bulk  cohesion 
energies  4.64  and  4.24  eV  for  the  diamond  and  FCC  structures,  respectively,  and  the  ab 
initio  potential  curve  of  the  silicon  dimer.  This  leads  to  the  constants  c\  =  0.225  eV. 
c’2  =  1.945  eV  and  c.3  =  —1.03  eV.  which  therefore  set  the  second  term  to  zero  for  Si?. 
Finally,  the  total  cohesion  energy  of  the  cluster  is  written  as  a  sum  of  the  band  structure 
energy  Ebs  and  the  repulsion  energy  Er, 


Cl 


#  -  t 


+  C3 


(4) 


Ecok  =  —(Ebs  +  Er) 


(5) 


which  is  maximized  to  find  the  structures  of  the  Si  clusters  in  the  present  work. 
The  dipole  transition  matrix  elements  (r)p?  can  be  written  as 

(r)p<7  —  }  '  Cgpa  Cp(ia  Rp  +  ^  {  C'7pj)Cppa(r)a6  • 
fia  ftab 


(0) 


where  the  atomic  transition  matrix  elements  are  (r)afc  =  JVj(r)  r  ^(r)  dr 


We  use  the  technique  of  the  one-electron  density  matrix  pP7  =  (a*ap).  which  ex¬ 
actly  takes  into  account  the  Fermi  statistics  of  electrons,  rather  than  the  conventional 
inanv-electron  matrix*6,2' .  Using  this  technique,  the  nlh-order  polarizability  is  con¬ 
ventionally  represented  by  means  of  double  Feynman  diagrams.  For  example,  one  of  the 
second-order  contributions  to  a12'  is  given  by 


where  Tpq  =  Tp  +  r?  is  the  pg-transition  linewidth.  Here  one  can  trace  the  diagram¬ 
matic  rules.  The  horizontal  lines  correspond  to  one-electron  states  over  which  sum¬ 
mation  is  implied.  These  lines  are  separated  by  vertices,  and  the  vertex  between  the 
lower  p  and  q  lines  corresponds  to  -e(r)7p,  and  between  the  upper  r  and  p  lines  to 
e(r)pr.  The  vertical  line  connecting  the  states  r  and  q  denotes  the  one-electron  propaga¬ 
tor  £>*  —  LJ7r  +  fT 7r]  .  where  W  is  the  sum  of  all  photon  frequencies  to  the  right 

of  the  propagator.  u:qr  =  (s7  —  £r)//i  is  the  transition  frequency  between  the  one-electron 

states,  and  the  free  p  lines  denote  pp0).  The  present  diagrammatic  technique  is  similar 
to  the  conventional  one,26,27  with  the  difference  here  that  all  the  operators  and  states 
are  one-electron  ones.  The  second-  and  third-order  polarizabilities  are  given  by  8  and  48 
diagrams,  respectively. 

NUMERICAL  RESULTS 

Maximization  of  Ecoh  (5)  yields  the  cluster  structures.  The  obtained  cluster  sym¬ 
metry  groups  and  important  electronic  structure  data  are  given  below  in  Table  1. 

TABLE  1.  Structure  data  for  the  Si?,  Si i o  and  Siu  clusters.  The  names  are  from  Ref.  17. 
Ecoh  (eV)  is  the  cohesion  energy  per  atom,  and  z.,  is  the  HOMO-LUMO  band  gap  (both 
in  eV  and  /an). 


Name 

y 

Point  exoiiD 

Ecoh 

( eV  /  atom ) 

-  7 

(eV, 

li  m) 

Coordination 

min.  max.  av. 

7 

3.S 

1.8 

0.69 

4 

6  4.6 

DBTA-I 

10 

C2v 

4.0 

1.4 

0.87 

4 

7 

TTP 

10 

3.9 

2.6 

0.4S 

3 

6  4.8 

TO 

10 

Td 

3.G 

2.9 

0.43 

3 

6  4.8 

13 

h 

4.4 

2.S 

0.44 

6 

12  6.7 

The  linear  polarizabilities  Re  a,;(u>)  (10-23  cm3'  of  the  clusters  obtained  for  dif¬ 
ferent  photon  energies  hu  are  given  in  Table  2.  It  follows  from  this  table  that  Rea  is 
nonsensitive  to  the  symmetry  of  the  cluster,  thus  bearing  little  structural  information.  In 
contrast,  the  cluster  absorption  spectra  (See  Fig.l),  governed  by  Ima.  strongly  depend 
on  the  cluster  structure:  the  lower  the  symmetry  of  the  cluster,  the  richer  its  spectrum 
with  more  red-shifted  offset. 

The  values  computed  for  the  second-order  polarizabilities  (10-28  esu)  = 

aS2\  — 2u\u,u)  and  0ijk{  —  u>,w)  =  a^2)(0;  —  ui,u),  governing  SHG  and  rectification,  are 
given  in  Table  3  (for  the  clusters  not  shown,  0  =  0  due  to  symmetry  requirements).  It  can 


section  (  1  O 


TABLE  2 


.V 
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ij  n~J  (eV): 
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1.77 
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2.33 
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5 

5 

0 

0 

15 
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G 
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i 

i 
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10 
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XX 

11 

12 

14 

16 
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3 

yy 

10 

11 

11 

12 

15 

3 

zz 

14 

14 

16 

IS 
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13 

10 

C3t) 
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12 

12 

13 

13 

14 

24 

z  z 

13 
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14 
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24 

10 

Td 

xx.yy,zz 

13 

13 

14 

14 
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20 

13 

h 

xx. yy, zz 

14 

14 

14 

15 

15 

17 

Photon  energy  (eV) 


Fig.  1.  Absorption  spectra  of  clusters  with  the  names  and  svmmetrv  groups  shown: 
(a)  Siio,  DBTA-I  (C2„);  (b)  Si10,  TTP  (C3„);  (c)  Si7,  ( Dsh );  (dj  Si10,  TO  ( Td ). 


TABLE  3 

N 

ijk  huj  (eV): 

o42*(  —  2u;:u;,u,’) 

0.0  10"2  0.5  1.0  1.17 

1.5 

a(2)(0;  —  u/',w) 

0.5  1.0  1.17 

10  C2v 

— 

-2.9 

-2.1  -3.9  5.3  3.9 

S.8 

-3.6  -8.0-13.8 

=yy 

-1.4 

-0.9  -1.6  1.4  -0.3 

2.1 

-1.7  -3.0  -4.2 

yyz,yzy 

-0.9 

-0.9  -1.5  S.l  3.5 

6.0 

-1.1  -2.0  -2.9 

zxx 

0.2 

0.3  0.8  -3.2  -5.1 

-2.4 

0.2  -0.0  -0.5 

xxz,  xzx 

0.4 

0.3  0.7 -3.1  6.2 

-45.6 

0.6  1.6  3.1 

10  Civ 

zzz 

l.S 

1.3  1.6  3.7  7.0 

31.2 

1.9  2.5  2.8 

yxx,xyx,xxy,  -yyy 

0.3 

0.3  0.3  0.6  1.0 

-4.7 

0.3  0.3  0.3 

zxx,zyy 

0.8 

0.3  0.3  0.4  0.5 

-16.0 

0.9  1.2  1.4 

xxz,yyz,xzx,yzy 

0.2 

0.3  0.4  0.9  1.7 

16.5 

0.2  0.3  0.3 

10  Td 

xyz ,  xzy, 
yxz , yzx. 
zxy,xyx 

-0.5 

-0.2  -0.3  -0.5  -0.8 

3.4 

-0.5  -0.6  -0.7 

be  seen  from  Table  3  that  for  the  lower  symmetry  clusters  Civ  and  C^v  the  magnitude  of 
3  is  considerably  higher  than  for  the  higher  symmetry  Td  cluster,  reaching  values  which 
are  characteristic  of  organic  molecules  with  high  nonlinear  responses. 

The  third-order  polarizability  7  =  o^3)  (data  are  not  shown)  is  symmetry-allowed 
in  all  the  cases.  Again,  the  maximum  responses  are  predicted  for  low-symmetry  clusters, 
reaching  high  values,  7  2  10~32  esu,  for  the  Civ  and  C73„  clusters. 

DISCUSSION 

Linear  polarizability  can  be  measured  for  single  clusters  in  experiments  using  laser 
light  scattering  and  photon-counting  detection.  However,  this  quantity  is  not  sensitive 
to  cluster  structure.  In  contrast,  the  optical  absorption  and  hyperpolarizabilities  do 
primarily  depend  upon  the  cluster  symmetry  and.  therefore,  bear  important  structural 
information. 

The  optical  absorption  of  single  clusters  in  jets  can  be  measured  similar  to  Ref.  30 
in  the  following  experiment.  Clusters  are  excited  by  tunable  probe  radiation  and  are 
also  subjected  to  powerful  IR  radiation  which  is  not  absorbed  from  the  ground  state  but 
ionizes  the  clusters  already  excited  by  the  probe  light.  The  clusters  ionized  are  detected 
by  the  mass  spectrometer. 

The  nonlinear  optical  responses  can  be  detected  if  it  is  possible  to  accumulate  clus¬ 
ters  in  host  media.  High  magnitude  of  the  polarizabilities  predicted  promises  possibility 
of  applications  in  optical  and  optoelectronic  devices. 
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